Two component Bose-Hubbard model with higher angular momentum states 
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Bose-Hubbard Hamiltonian of cold two component Bose gas of spinor Chromium atoms is studied. 
Dipolar interactions of magnetic moments while tuned resonantly by ultralow magnetic field can 
lead to a transfer of atoms from the ground to excited Wannier states with a non vanishing angular 
orbital momentum. Hence we propose the way of creating of P x + iP y orbital superfluid. The spin 
introduces an additional degree of control and leads to a variety of different stable phases of the 
system. The Mott insulator of atoms in a superposition of the ground and vortex Wannier states 
as well as a superposition of the Mott insulator with orbital superfluid are predicted. 



I. INTRODUCTION 

Ultracold atoms provide a playground for mimicking 
condensed matter and studying novel quantum many- 
body phenomena 0, 0. Recently, there has been par- 
ticularly impressive progress in two areas of physics of 
ultracold atoms: the area of ultracold dipolar gases 
and the physics of orbital lattices 0- In this paper we 
combine these two areas and explore the effect of two- 
body dipolar interactions of magnetic atomic moments 
in a lattice potential. We study dipolar gases in their full 
complexity including spin as a dynamical variable (as op- 
posed to be a conserved quantity), and magnetic dipolar 
interactions coupling different orbital states of involved 
magnetic components. That introduces additional phys- 
ical processes into play and new degrees of control to the 
standard Bose-Hubbard model. 

Spinor gases in a lattice have been studied in the con- 
text of Mott insulator (MI)- superfluid (SF) transition 
Jj|. In general dipolar interactions lead directly to the 
dynamics of spin degree of freedom but up till now in 
lattice systems this phenomenon was neglected, i.e. it 
was assumed that spin is frozen Q. In such situations 
electric and magnetic dipoles are practically equivalent 
- they introduce long range correlations. On the other 
hand it is known from studies of gases confined in har- 
monic traps in the mean field limit 0, H| that taking the 
dynamics of spin into account may modify properties of 
the ground state of the system. 

Spin dynamics may result from contact or dipolar in- 
teractions. In the former case the total spin of inter- 
acting atoms remains unchanged (magnetization of the 
sample is constant). Qualitatively different phenomena 
take place when spin dynamics is triggered by the dipo- 
lar forces. The atomic magnetic moment originates from 
the spin which contributes to the total angular momen- 
tum of the system. When magnetic dipole changes due 
to the dipolar interactions, its variation must be accom- 
panied by corresponding dynamics of the orbital angular 
momentum. Magnetic interactions can lead to a trans- 



fer of angular momentum from spin to orbital degrees 
of freedom. This phenomenon, discovered in ferromag- 
netic solid samples, is known as Einstcin-dc Haas effect 
Not a long range character of magnetic dipolar 
interactions but rather their relation to the angular mo- 
mentum plays a crucial role in this phenomenon. This 
makes a fundamental difference between magnetic and 
electric dipoles. 

The main issue of our study is to account for the spin 
degree of freedom in the lattice environment. Spin flip- 
ping processes in the lattice could lead to an appearance 
of the orbital P x + iP y superfluid. Recently orbital su- 
perfluids were created in experiment The authors 

utilized a resonant tunneling in a particularly designed 
lattice potential. 

In this paper we show another way of creating orbital 
superfluid by means of the resonant Einstein-de Haas ef- 
fect. The atom which flips its spin has to gain some ad- 
ditional kinetic energy necessary to support its rotation. 
This energy is typically much larger than the energy of 
dipolar interactions and conservation of energy strongly 
suppresses the spin dynamics. The transfer of atoms be- 
tween two spinor components can be enhanced by tuning 
energies of states involved via Zeeman effect [rlfl . We 
extend this idea to lattice gases. Dipolar effects signif- 
icantly modify the MI-SF transition lead to new phases 
of the system with quantized vortices in MI or/and SF 
regimes. 

The paper is oranized as follows: in Section II we in- 
troduce the two component Bose-Hubbard model with 
dipolar interactions coupling different Wannier states, in 
Section III we present a phase diagram for the system 
while in Section IV we discusse validity and limiatations 
of the model. 



II. THE MODEL 

We assume that Cr atoms are in a 2D optical square 
lattice. To fix the parameters we consider a realistic sit- 



2 



uation of the lattice described by the periodic potential 
Vb [sin 2 (2ttx) + sin 2 (27ry)] Here A = 523 nm is the wave- 
length of light beams creating the lattice and Vq is the 
barrier height. A characteristic energy of the problem, 
i.e. the recoil energy is E r = h 2 (2tt) 2 / (2mA 2 ) . We ex- 
press all energies and lengths in units of E r and A respec- 
tively. Confinement along the z direction is provided by 
a harmonic potential mui 2 z 2 /2 of frequency huj z = 16E r . 
At each lattice site we choose two wave functions cen- 
tered at the given site (xi,yi) to form a single particle 
basis of the two component system. The basis allows 
to account for the resonant transfer of atoms between 
ras = 3, I = and 7715 = 2 and I = 1 states in the pres- 
ence of magnetic field aligned along the z-axis. The low- 
est energy state ip a (x, y, z) ~ W (x)W {y) exp(-z 2 oj z /2) 
is effectively coupled to the excited state with one 
quantum of orbital angular momentum ipb(x,y,z) ~ 
[Wi(x)Wo(2/) + iWo(a;)Wi(2/)]cxp(-z 2 w z /2). The state 
is a single site analogue of a hamonic oscillator state 
- (x + iy)exp[-{x 2 +y 2 )/2-z 2 uj z /2}. W (x) andWi(.T) 
are the ground and the first excited Wannier states in 
a ID periodic potential of the form Vq sin 2 (27rx). Single 
particle energies of the two essential states are denoted 
by E a and Eb respectively. 



Limiting the subspace of essential states is a crucial 
approximation in our study. It is possible only due to a 
weakness of dipolar interactions. In fact there are several 
channels of binary dipolar collisions leading to different 
excited Wannier states. However, we can choose the de- 
sired channel by a proper adjustment of the resonant ex- 
ternal magnetic field . Typically the energy difference 
between atoms in the ground and in the excited Wannier 
states is much larger then dipolar energy which is the 
smallest energy scale in the problem (except vanishing 
tunnelings case), Ed ip = 10~ 4 E r <C Eb — E a ~ E r . How- 
ever, at resonant magnetic field Bq, E a — gfJ-sBa = Eb, 
the two energies are equal and the spin transfer between 
the components becomes efficient on a typical time scale 
fr/Edip — 10 _2 s. Here /is is the Bohr magneton and 
g = 2 is the Lande factor. Only then the system can 
dynamically redistribute particles between the two com- 
ponents without violating energy conservation. A char- 
acteristic width of the resonances is small [lH, of the 
order of Edi P ~ gfisB, i.e. B « 100/iG. We assume that 
no other states can be effectively coupled (see a more de- 
tailed discussion of the validity of this model in Section 
IV). 



In effect a two-component system is realized with a- 
componcnt corresponding to atoms in 7715 = 3 and I = 
state while atoms in 6-component have mj = 2, I = 1. 
Single site basis states are \n a , rib), where n c is a number 
of atoms in c-componcnt (c = a, 6). The Hamiltonian of 



the system is: 
H = [( E a - g^BB) a\a t + E b b\h + U ab a\b\a t b t 
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The parameters depend only on lattice height Vq and 
confining frequency ui z in the z-direction. U a , Ub, U a b are 
the contact interaction energies plus the part of dipo- 
lar energy which has the same form as corresponding 
contact term, D is the on-site dipolar coupling of the 
two components, while J a and Jb are tunneling energies. 
The Hamiltonian (p} is an interesting modification of the 
standard Bose-Hubbard model. 

The on-site contact interactions U a ,Ub, and U a b cannot 
change a total spin [l6j . Dipolar two body interac- 
tions are much smaller than the contact ones; we keep 
only those dipolar terms which lead to a spin dynam- 
ics. Moreover, only on-site dipolar effects are accounted 
for in the Hamiltonian ([1]). Dipolar potential, although 
long range, is so weak that we can ignore dipole-dipole 
interactions between atoms at neighboring sites in the 
considered range of small tunnelings. 

Unlike tunneling between ground Wannier states J a , 
the tunneling energy Jb of the excited state is negative 
because the wave function the ijjb(x, y, z) is antisymmet- 
ric in x and y. Therefore the state with 'antifferomag- 
netic' order of phases between neighboring sites has lower 
energy than the state where phases of the exited Wannier 
functions are the same. For the opposite on-site phases 
of the excited Wannier states both J a and Jb are positive. 
This case is considered here. 



III. PHASE DIAGRAM OF THE MODEL 

We limit our study to a small occupation of a lattice 
site: not more than one particle per single site on average. 
The resonant magnetic fields equilibrates single particle 
energies of states |1,0) and |0, 1), i.e. Eb = E a — g^sBo- 
E a and Eb depend on the lattice height thus the resonant 
magnetic field varies with Vo, B = B (V~o). 

Even with a single particle per site the dipolar inter- 
actions couple ground and excited Wannier states due 
to the tunneling in a higher order process. The trans- 
fer between |1,0) and |0,1) states is a sequence of: 
adding an atom to the a-component at a given single 
site |1,0) — > |2,0) via tunneling, followed by the dipolar 
transfer of both a-species atoms to the excited Wannier 
state |2,0) — > |0, 2), and finally the tunneling which re- 
moves one 6-component atom from the site |0, 2) — > |0, 1). 
The two considered states are therefore coupled provided 
that tunneling is nonzero. 

Now, following the standard mean field approach of 
Fisher et al. 17] we find thcrmodynamically stable phases 
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FIG. 1. Phase diagram for 2D square lattice at the resoance, 
2 = 4. The regions are: M - Mott insulator with one parti- 
cle in equal superposition od a and fa states, MS - superfluid 
in a and fa components (fa-dominated) and Mott insulator in 
the orthogonal supperposition, S - superfluid phase of su- 
perposition of a and b components, Sb - superfluid in the 
fa-component. In the inset the diagram for 2 = 3 together 
with chemical potential fl(N) for a given number of particles 
obtained from the exact diagonalization. The lines, from bot- 
tom to top correspond to occupation equal to TV = 2, . . . , 9 
as indicated. For /i > Ub (light grey region) the ground state 
of the system is a two particle state, therefore in this regime, 
the phases shown are thermodynamically unstable. They are 
stable, however, with respect to one particle hopping. 



of the system in the choosen subspace. The Hamilto- 
nian (p} is translationally invariant, we assume the same 
property is enjoyed by the lowest energy state. Intro- 
ducing superfluid order parameters for both components: 
4>(a) — ( a i) an d 0(6) = (pi) as well as the chemical po- 
tential [i, the Hamiltonian of the system can be approx- 
imated by a sum of single site Hamiltonians Hq + Hj 

H = -/i(a f a + fa f 6) + ^U a o)a)aa + ^Ubtftfbb 

+U ab a%*ab + D(tftfaa + a*<Jbb), (2) 
Hi = -zJ a <j>*a)a- zJb<p*qb + h.c. (3) 

Notice we skipped indices enumerating sites. In ([3]) z is a 
number of neighbors and depends on the lattice geome- 
try. For a 2D square lattice z = 4. Hamiltonian H + Hi 
does not conserve number of particles: it describes a sin- 
gle site coupled to a particle reservoir. Order parame- 
ters <f>r a \ and vanish in the MI phase and hopping 
of atoms is suppressed. Only in the SF regime number 
of particles per site can fluctuate. Close to the bound- 



ary, on the SF side, 



and 



(6) 



can be treated as small 



parameters of the perturbation theory. 

The single site ground state becomes unstable if the 
arc different than zero. The 



mean field 



(6) 



consistency condition: 
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where c = a,b. In the lowest order of the perturbation in 
the order parameters, the set of equations ([4]) becomes 
linear and homogeneous. Vanishing of its determinant is 
a necessary condition for nonzero solutions for rf>r c \ . This 
condition determines lobs shown in Fig. Q] 

In the low temperature limit (/3 — > oo) the partition 
function reduces to a single lowest energy state contri- 
bution Z(/3) = e- pEo . The energy Eq depends on the 
chemical potential fi. Moreover, for fj, < Ub < U a the 
only contribution to Eq. (0} comes from eigenstates of the 
Hamiltonian with zero, one and two particles. Our anal- 
ysis is Unfitted to this case only. 

For negative chemical potential, \i < 0, the single site 
ground state is |0,0) vacuum state (dark grey region in 
Fig- [J). With increasing tunneling (and fixed fj,) parti- 
cles appear in the superfluid vortex 6-phase (Sb)- Only 
at larger tunnclings some atoms do appear in the a- 
componcnt and both: 'standard' and P a 
supcrfluids coexist (S). 

Situation becomes more complicated for larger chemi- 
cal potential < fj, < Ub- At the resonance, B = Bo, the 
ground state is degenerate if tunneling is neglected: the 
states 1 1 , 0) and |0, 1) have the same energy, Eq = —fj,. 
The degeneracy is lifted via tunneling in the second or- 
der of the perturbation. In addition a position of the 
resonance is shifted towards smaller magnetic field val- 
ues. Analysis of the effective Hamiltonian (compare 
[lSl] ) indicates that in the resonant region the single 
site ground state is a supperposition of both compo- 
nents \g) = ai|l,0) — a^O, 1). Exactly at resonance 
ot\ = q-2 = 1/a/2- While crossing the resonance the 
ground state switches from |1,0) to |0,1). The width 
of the resonance AB can be estimated perturbatively 
to be gpL B \AB\ w lO" 6 ^ for V = 25E r while for 
lower barriers, Vq = l0E r , the resonant region is broader 
5^b|AB| w lO -3 ^. Due to its small width the reso- 
nance can be hardly accesible particularly for small tun- 
nclings. Away from the resonance the standard phase 
diagrams for a or 6 component emerge. 

In Fig. [T] we show regions of stability of different pos- 
sible phases of the system at resonance i.e when \g) = 
(|1, 0) — |0, 1))/ y/2. For small tunnelings the system is in 
the Mott insulating phase (M) with one atom per site. 
Every atom is in the superposition of the ground and 
the vortex Wannier state. At the blue line, the border 
of (M) lobe, Eqs.(f4]) allow for nonzero solutions for <j)( a \ 
and 4>fh\. Eqs.(f4]) become diagonal if Hi is expressed 
in terms of bosonic operators — (n a a) + Kbb') and 



flt = (. 



KbCl' 



K a b') where 



1 and both cocfi- 



mean fields can be obtained numerically from the self- 



tients of the superposition depend on the tunellings J a 
and Jb- The operators create an atom in two orthogo- 
nal superpositions of a and b states. At the border of 
the Mott phase (M) the mean value of the operator B is 
different from zero and a nonvanishing superfluid compo- 
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nent, ^ b = —i^b4>(a)+ K -a4>(b), appears in the (MS) region. 
Our numerical results show that K a ~ —0.99 and the ra- 
tio {Kb/no) 2 — 0.02 is small at the edge of stability of the 
Mott insulator. Therefore ~ b\ i.e. the supcrfluid 
b is dominated by the orbital 6-component. The mean 
field corresponding to the ~ operator is zero in the 
discussed region. The system is therefore in equal su- 
perposition of the Mott insulating and supcrfluid phases. 
The Mott phase is dominated by the a-component and 
the superfLuid phase is overwhelmed with the 6-species. 
Both components, however, contain a small minority of 
remaining species. 

At larger tunneling the system undergoes another 
phase transition as Eqs.Q allow for another nonzero 
mean field. Now the mean value of A departs form zero 
defining the border of the 'bigger' lob. Mott component 
of the ground Wannier state becomes unstable. The ad- 
ditional mean field \1/ a = ^a0(a) + K &0(fi) appears in the 
(S) region. Again n a — 0.97 and the maximal value of 
{Kb/K a ) 2 ~ 0.06 is small. The a-species dominate the 
if? A superfluid component. Both if? a and if?B superfluids 
exist in the (S) region. 

All the above findings are supported by direct inspec- 
tion of the true many body ground state obtained by 
exact digonalization of the many body Hamiltonian in a 
small 2x4 rectangular plaquctte with periodic boundary 
conditions for total number of particles N = 1, . . . , 10. 
Note that each site has three neighbors, z = 3, in this 
case. Resonance condition is reached by finding the mag- 
netic field for which both a and 6 species are equaly popu- 
lated. Calculations for z = 4 require much larger number 
of sites and are numerically unreachable. In the inset of 
Fig. [T] we compare the exact results with the mean field 
ones but for z = 3. The lines in the inset correspond to 
the constant number of particles per site obtained from 
the relation n(N) = [E {N + I) - E a (N - 1)] /2. They 
allow to trace the phases the system enters while adia- 
batically changing the tunneling at fixed particle number. 
The (M) and (MS) phases can be reached with one par- 
ticle per site only (8 particles in the plaquette). Direct 
inspection of a structure of the many body ground state 
fully confirms the stable phases of the system described 
above. In particular the ground state in the (MS) region 
can be approximated (with the accuracy of about 4%) by 

tsM - wherc I") is the vaccum 

state. 

In addition we calculated a hopping, i.e. the mean val- 
ues of the following hopping operators : h a = Yl/j) ( a ] a i) 
and hi, = (bjbi)- These operators annihilate a par- 
ticle at a given site and put it in a neighboring site. 
They might be viewed as number conserving analogons 
of the mean fields 0( a ) and 4>{b)- In Fig. [2] we show 
the hopping for the case of one particle per site. For 
large tunnelings both a and b hopping are large - the 
components arc in the supcrfluid phase. Entering the 
MS phase, J a /U a ~ 0.064, the hopping of a-component 
rapidly falls down while hopping of 6-atoms remains big 




FIG. 2. Hopping for the lowest energy state in a 2 x 4 pla- 
quette obtained from the exact diagonalization. Upper line - 
b component, lower line - a component. 



- the system enters a-component dominated Mott insu- 
lator superimposed with 6-component dominated super- 
fluid. At J a /U a — 0.002 both hoppings tend to zero - the 
system enters the Mott phase with equal occupation of 
both species. This confirms results based on the Fisher 
method. 



IV. VALIDITY OF THE MODEL 

Finally let us discuss possible limitations of the valid- 
ity of the model discussed above. As we study a sta- 
bility of the Mott phase we consider the case of deep 
optical lattices where tunneling is a small perturbation 
only. It is very natural to assume that dipolar interac- 
tions couple the ground Wannier state to the orbital state 
at each lattice site, and the system posseses the transla- 
tional symmetry. Moreover, we have assumed that locally 
the potential at a given site has almost perfect axial sym- 
metry with respect to the site center. Therefore, the local 
site Hamiltonian preservs projection of the total angular 
momentum, and the only state coupled to the ground 
Wannier one is of the type ~ (x + iy), where x and y 
are measured with respect to the site center. This state 
is the eigenstate of the projection of the orbital angular 
momentum on the z-axis. 

Three comments arc in order. 



A. Role of anharmonicity. 

Due to high selectivity of magnetic resonances we have 
a freedom of choosing a given channel of dipolar colli- 
sion by a proper adjustment of the external magnetic 
field. In particular, we study the channel where the z- 
comonent of the relative orbital angular momentum of 
interacting particles changes by two quanta, A.L Z = 2. 
Assuming that each of two colliding atoms are initially 
in the spherically symmetric ground state, the lowest en- 
ergy final state of the two atoms has a form | vortex) ~ 
(xi+iyi) 2 + {x 2 +iy2) 2 -'2{xi+iyi){x2+iy2)- Note that in 
the harmonic trap of radial frequency w the state | vortex) 
corresponds to a superposition of two states: |v2) ~ 
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(xi + iyi) 2 + (x 2 + ij/2) 2 and |vi) ~ (xi + iyi)(x 2 + iyz), 
where Xi, yi are particles coordinates. For |v2) one of 
colliding atoms aquires two quanta of rotation while the 
second atom remains in the spatial ground state, i.e. the 
energy of |v2) is E2 = 2hu. On the other hand |vi) 
represents the situation where each of two atoms gets 
one quantum of rotation resulting in the total energy 
2E\ = 2fko. Evidently for equally spaced harmonic en- 
ergy levels both states are degenerate, E 2 = 2E\ and 
both the conservation of angular momentum and the con- 
servation of energy can be satisfied. 

The situation becomes different in the optical lat- 
tice because of anharmonicity of the lattice potential. 
The state |v2) has energy of the second Wannier state, 
Ei- This energy is smaller than twice the energy of 
the first excited Wannier 2E\ of the state |vi). Even 
for high barriers, i.e. Vq = A0E r the energy splitting 
2Ei — E 2 = 10 _1 -E r is significantly larger then dipolar en- 
ergy of Cr atoms, Edi P = 10~~ 4 E r . Therefore, by means of 
magnetic field tuning one may select the resonant trans- 
fer of atoms due to dipolar interactions bringing both 
interacting particles to the state with one quantum of ro- 
tation while making the transfer to the second Wannier 
state nonresonant (and not efficient). This is a situation 
considered in the present paper. Note that a proper ad- 
justment of the magnetic field may make the excitation 
of D orbital |v 2 ) resonant - the situation not considered 
here. Anharmonicity of the lattice potential, although 
small, plays thus an important role. 



is not strictly true. The single site potential cannot be 
approximated by a harmonic one if fine details are to 
be studied. In a sqaure 2D lattice every site has four 
neighbours and the square symmetry of the lattice in- 
fluences the single site potential. The quartic terms in 
the expansion of Vo(sin 2 (27ra;) + sin 2 (27ry)) potential are 
relevant. For this reason the two particle state, |vortex) 
corresponding to Lz = 2 is not the eigenenergy state of 
single particle plus contact interaction on-site Hamilto- 
nian. It is a superposition of three two-particles states 
of different energy instead. The fine structure results 
both from the anharmonicity and the anisotropy of the 
trapping potential. The anisotropy of the trap cannot 
be reduced even for very high lattices. We checked that 
even for Vq = 100E r the energy splitting is significantly 
larger then dipolar interaction energy. Large magnetic 
moments, leading to larger dipolar energy could help to 
overcome this problem. Conservation of energy allows to 
tune independently only to the one of the three compo- 
nents of the jvortex) state. Weak dipolar interactions re- 
solve this fine structure of two-particle energy states. To 
observe the Einstein de Haas effect in optical lattices one 
should use the lattice geometry for which the anisotropy 
due to the lattice symmetry is substantially reduced. To 
this end a 2D triangular lattice with every site having 6 
neighbours might be promissing. The other way out is 
to rotate every lattice site around its axis similarly as in 
the experiment [l9j |. 



B. Vorticity versus tunneling 

The second issue is related to the tunneling of vortex- 
like states. As the lattice states have C4 symmetry the 
angular momentum need not be conserved in the tunnel- 
ing. With our choice of alternating phases of the excited 
Wannier states the tunneling coefficient in the excited 
band is positive, Jb > 0. Therefore tunneling of the right 
handed vortex ~ (x + iy) to the vortex of the same vor- 
ticity at neighbouring site is equal to i_r? = Jb + J a and is 
larger then tunneling with simulataneous change of the 
vorticity, i.e. to ~ (x — iy) state, ti = Jb — Ja- The 
difference is small since tunneling in the lowest Wannier 
state J a <C Jb but significant. The tunneling decreases 
the system energy thus the larger tunneling for the pro- 
cess preserving vorticity will decrease the system energy 
more and will be preferred. This observation allows for 
including only right handed vortex in our single particle 
basis and omitting the left handed one. 



C. Single site anisotropy 

In our model we have assumed that the single site po- 
tential is isotropic, i.e. the two particle state produced 
in the dipolar interactions has both the well defined en- 
ergy and the relative angular momentum. In fact this 



V. CONCLUSIONS AND OUTLOOK 

In this paper we studied the model Bose-Hubbard sys- 
tem with two Wannier states in optical lattice. We show 
that weak dipolar interactions can be resonantly tuned 
to couple the ground Wannier state to the excited one 
with higher orbital angular momentum. We have stud- 
ied a case of at most one particle per site on average. 
Even in this case, we predict various novel phases of the 
system. The phase diagram of the system significantly 
depends on the magnetic field. On the resonance we pre- 
dict three distinct phases of the system: i) the Mott in- 
sulator of superposition of ground and vortex states ii) 
the a-component dominated Mott superimposed with b- 
componcnt dominated superfluid, iii) two superfluids in 
particular combination of both species. We also discuss 
some limitations of our approach stressing that harmonic 
approximation has to be used with caution when study- 
ing orbital physics in optical lattices. 

Higher densities (more particles per site) are more fa- 
vorable for dipolar transfer, the related physics will be 
discussed elsewhere. It is worth noting that our results 
may be direcly related to the very recent experiments, 
in which spin relaxation in an ultracold dipolar gas in 
an optical lattice was observed in a presence of ultra low 
magnetic field [l|l[2(|. Although, so far, no vortices have 
been found, we hope that the present work will help to 
identify the regime of parameters, in which generation 
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of P x + iPy supcrfluid and appearance of novel quantum 
phases occurs. 
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